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Abstract 

By a classical theorem transversal homoclinic points of maps lead to 
shift dynamics on a maximal invariant set, also referred to as a homoclinic 
tangle. In this paper we study the fate of homoclinic tangles in parameter- 
ized systems from the viewpoint of numerical continuation and bifurcation 
theory. The bifurcation result shows that the maximal invariant set near a 
homoclinic tangency, where two homoclinic tangles collide, can be charac- 
terized by a system of bifurcation equations that is indexed by a symbolic 
sequence. For the Henon family we investigate in detail the bifurcation 
structure of multi-humped orbits originating from several tangencies. The 
homoclinic network found by numerical continuation is explained by com- 
bining our bifurcation result with graph-theoretical arguments. 
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1 Introduction 

We consider parameter dependent, discrete time dynamical systems of the form 

Xn+l = f{Xn,X), G (1) 

where /(-, A), A G R are smooth diffeomorphisms in R'^. We assume that the sys- 
tem ([1]) has a smooth branch of hyperbolic fixed points and our main interest is 
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in branches of homoclinic orbits that return to these fixed points. Generically one 
finds turning points on these branches which correspond to homochnic tangencies 
where stable and unstable manifolds of the fixed point intersect nontransversally, 
see [15], p] for the precise relation. While the dynamics near transversal in- 
tersections are well understood through the celebrated Smale-Shilnikov-Birkhoff 
Theorem (see [26], [25], [9], [21]), the picture near homoclinic tangencies still seems 
to be far from being complete. Among the many references, we mention the 
monograph [20] which contains a detailed geometrical study of the bifurcations 
that occur near tangencies, the work [7] which supports the generic occurrence 
of homoclinic tangencies of all orders, and the paper [17] which proves shift dy- 
namics in arbitrarily small neighborhoods of the tangency. We also mention 
that homoclinic orbits with a tangency can be computed numerically in a robust 
way by solving boundary value problems on a finite interval, and that the errors 
caused by this approximation have been completely analyzed, see [TBI [T5| 13]. 

In this paper, we consider bifurcations of tangential homoclinic orbits from a 
local as well as from a global viewpoint. 

The local study determines the elements of the maximal invariant set in a 
neighborhood of the tangent orbit and of the critical parameter from a set of 
bifurcation equations. Using the same shift space as in the transversal case, 
we associate with any sequence of symbols a bifurcation equation that describes 
those branches of orbits which have the return pattern of the symbolic sequence. 
Such a result does not fully resolve the dynamics near tangencies, but reduces the 
problem to a set of perturbed bifurcation equations for which the unperturbed 
form is known (similar to Liapunov-Schmidt reduction). For example, multi- 
humped homoclinic orbits that enter and leave a neighborhood of the fixed point 
several times, relate to a perturbed system of hilltop bifurcations, cf. [0]. Our 
main results will be stated in Section |2] with the proofs deferred to Sections O |6l 

The global viewpoint asks for possible bifurcations of multi-humped orbits 
that are known to emerge from the tangencies. We take the Henon family as a 
model equation for a detailed numerical study of the homoclinic network that 
arises from a total of 4 primary homoclinic tangencies. It turns out that the 
connected components of this network are by no means arbitrary. Rather, they 
follow certain rules governing the bifurcations of multi-humped orbits. Combining 
these rules with graph-theoretical and combinatorial arguments allows to predict 
the structure to a large extent. Only some fine details are left to numerical 
computations as will be demonstrated in Sections [3] and ID 

2 Setting of the problem and main results 

The aim of this section is to state our main result on bifurcation equations near 
homoclinic tangencies. We first describe the setting and state our assumptions: 

Al / e C°°(R'' X Ao, BI') for some open set Aq C R and /(-, A) is a diffeomor- 
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phism for all A G Aq, 

A2 /(^(A), A) = ^(A) for some smooth branch ^(A) G R^ A G Aq, 

A3 /x.(^(A), A) G R''''' is hyperbolic for all A G Aq. 

Clearly, if ,^0 is a hyperbolic fixed point of /(■, Aq) for some Aq G R then A2 and 
A3 follow for some neighborhood Aq of Aq. Replacing / by g{x, A) = /(x + ,^(A + 
Aq), A + Aq) — ^(A + Aq) shows that A2, A3 can be assumed to hold for the trivial 
branch ,^(A) = and for a neighborhood Aq of zero. This will be our standing 
assumption throughout Sections [5] and [61 

It is well known that transversal homoclinic orbits lead to chaotic dynamics 
on a nearby invariant set commonly referred to as a homoclinic tangle. Let us 
first assume that this situation occurs at some parameter value A G Aq. 

A4 For some A G Ag there exists a nontrivial homoclinic orbit — (2^n)n6S; 
i.e. lim^^iooa^n = ^ '■= ^(A) and x„ 7^ ^ for some n G This orbit is 
transversal in the sense that the variational equation 

Un+i = fx{xn, ~X)yn, n G Z (2) 

has no nontrivial bounded solution on 

In this case the stable and the unstable manifold of ^ intersect transversally 
at each Xn and the set H = {xn}ne'z U is hyperbolic, cf. [22]. Moreover, there 
exists an open neighborhood U of H such that the dynamics on the maximal 
invariant set 

M{U,~X) = {x eU : r{x,~X) eWineZ} (3) 

is conjugate to a subshift of finite type (see the Smale-Shilnikov-Birkhoff Homo- 
clinic Theorem in [9] and [2T1 Chapter 5] for a proof). To be precise, let iV > 2 
and let 

5^ = {0,l,...,iV-lf 
be the shift space with symbols which is compact w.r.t. the metric 

d{s,t) = ^2-^%sj -tj\, s = {sj)j^^ E Sn. (4) 

Let /3 be the Bernoulli shift 

P{s)i = Sj+i, i G s e Sn- 
Consider a special subshift of finite type, see [18] 

n^ = {seSN: Ai^X^^ = 1 G 
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generated by the N x N binary matrix 

/l 1 
1 



V ••• 





1 

0/ 



G{0,1} 



NxN 



Then there exists a neighborhood U of H, an integer N > 2 and a homeomor- 
phism h : Qn ^ M{U, A) such that 



f{-,~X)oh = ho (3 inn 



N- 



(5) 



A continuation of the transversal homochnic orbit w.r.t. the parameter A 
leads to a curve of homoclinic orbits that typically exhibits turning points. As 
an example we refer to Figured] for the Henon map. Parts of the branch that can 
be parametrized by A belong to transversal homoclinic orbits while (quadratic) 
turning points correspond to homoclinic orbits with a (quadratic) tangency, see 
Theorem |3] for a precise statement. In this case we replace Assumption A4 by 

B4 For some A G Aq there exists a nontrivial homoclinic orbit — (2^n)nes 
converging towards ^ = ^(A). The orbit is tangential in the sense that the 
variational equation 



Vn+l = fx{Xn, X)yn, U ^TL 



(6) 



has a non-trivial solution u% = [unju^z that is unique up to constant mul- 
tiples. 

Since the fixed point stays hyperbolic we have exponential decay for both the 
orbit and the solution of 



i.e. for some a, Cg > 



\u 



nil ^ C*ee 



Therefore, we may normalize 



T 



(7) 



(8) 



In the following we use (■, ■) to denote the inner product in i^. The assumption 
on ([6]) in B4 holds if and only if the tangent spaces of the stable and unstable 
manifold have a one-dimensional intersection, i.e. 

T,,^W'{0 n T,„iy"(e) = span(M„), n e 
We refer to Theorem [3] and to [151 Appendix] for a more general statement. 
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Consider open neighborhoods U C IR!^ of H = {xn}ne'z U {^} and A C Aq of 
A, respectively. Our main interest is in the dynamics on the maximal invariant 
set 

M{U, A) = {{x,X) eU X A: /"(a;, X) E U e TL}. 

As in the transversal case, we will still work with the subshift (I^at,/?) but the 
conjugacy ([5]) will be replaced by a set of bifurcation equations. For any s 
define the index set 

/(s) = G Z : s„ = 1}, (9) 
and note that / : f^Tv — ^ Z(iV) C 2^ is bijective, where 

Z(iV) = {J -k\>N Vj, keJ,3^ k}. (10) 

With any s G f27v we associate the Banach space 

£^{s) = {r G R^('^ : ||r||oo < oo}, 

llrlloo = sup |r(^)|, 5, = {r G : ||r||oo < p}- 

Our aim is to determine the elements of M{U, A) from a set of bifurcation 
equations 

gs{T,\) = 0, reB,^,\eA, (11) 
where s E Qn, Pt > is independent of s and 

(r,A) ^ gs{rA) 

is a sufficiently smooth map. Note that ( |TT|) constitutes a finite or an infinite 
system of equations depending on the cardinality of /(s). 

In order to formulate the precise statement we define the pseudo orbits 

l&I{s) 

Equation ([7]) shows that ^^(s) is a bounded sequence, in particular 

\\Pn{s) - el <C = Cj-^^, nETL. (13) 

1 — e ° 

Setting = ^ for all n G Z we write ( IT2l) more formally as 

Here and it what follows we use the symbol /3 to denote the shift of sequences in 
R*^. Thus /3 acts as an operator in sequence spaces such as &'{^\ 1 < p < oo. 
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Similarly, for every r G i°°{s) we define the bounded sequence 

vzis,T) = (14) 

iei{s) 

Note that the sequence P'zis) has humps at the positions defined by /(s) and that 
P'z{s) is a pseudo orbit of f{-,X) with a small error, see Lemma W2\ The term 
V'z{s, t) shifts the solution of the variational equation to the positions defined by 
J(s) and combines them linearly. 

Theorem 1 Let assumptions Al - A3 and B4 hold. Then there exist constants 
< Tr < Pt, N & ¥i and neighborhoods U of H, A of X and for any s E Qn 
smooth functions 

x^,s : Bp^xA^ i'^iR^) 

with the following properties. 

(i) For any point {yo, \) G M{U,A) with orbit yn = f"'{yo,X), n E 1i there 
exists an index v eTL and elements s G ^at, r G Bp^ C £°°(s) such that 

(^''yw. = Xz,s{t^ ^) + Ps(s) + vz{s, r), (15) 
9sir,X) = 0. (16) 

(a) Conversely, if s G Qn, t G Br^ C i°°{s), X G A satisfy f|T6|) . then there 
exists u E 1i such that {yn, A)„ga. with given by f|T5|) . belongs to M{U, A). 

Remark 2 Theorem U\ reduces the study of M{U, A) to the set of bifurcation 
equations (fT6|l with a symbolic index s E Qn- It may be regarded as a type of 
Liapunov- Schmidt reduction though we have not formally put it into this frame- 
work. The construction of the neighborhood U x A uses some features from the 
transversal case f2l[ Theorem 5.1], but is considerably more involved, see Sections 
andl^ We also note that we were not able to prove that one can take r-r = 
which would give a complete characterization of M{U, A) in terms of f|T5|) . f|T6l) . 
Another issue which has not yet been resolved, is continuous dependence of the 
functions Xs,s (^nd gs on the symbolic sequence s with respect to the metric (jlj). 

The functions gs and Xs,^ have several properties that we discuss next. 
Due to B4 the adjoint equation 

i/r+i/x(xn+i, >) = yl, nez (17) 

has a non-trivial solution u;^ that is unique up to constant multiples, cf. [221 
Section 2]. It decays exponentially as in ([7]) and can thus be normalized such 
that I Ill's 1 1^2 = 1. Without loss of generality we take Cg in ([7]) such that 

< Cee-"l"l, nGZ. (18) 
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As is shown in [15] the quantities 

Ca = {W^, {fx{Xn, A))„ea), Cx = ^{W^, {fxxiXn, A)u^)nes) (19) 

characterize the behavior of the branch of homochnic orbits which passes through 
(^z, A). 

Theorem 3 The operator F : £°°(R^') x R ^ £°°{R'') defined by 

F{x^, A) = {{xn+i - f{xn, A))„ez) (20) 

has a limit point at (xg, A) in the sense that -F(xs, A) = and 

Af{D^F{xz,)\)) = span{Ms}. 

The limit point is transversal, i.e. 

DxF{x^, A) ^ n{D,F{xz, A)) if and only if cx + 0. 

Moreover, it is a quadratic turning point, i.e. 

DlF{xz, A)n| ^ n{D,F{x^, A)) if and only if c,. ^ 0. 

Remark 4 For transversal homoclinic orbits (X ^ X), the Sacker-Sell spectrum, 
cf. l2^ , of the variational equation ^ is a pure point spectrum 

SED(A) = {|/i|:/iGa(/,.(e(A),A))}. 

At a turning point, we find a spectral explosion to a continuous Sacker-Sell spec- 
trum 

Sed(A) = : ji e (T{fxi^,X))} U 

where 

Us = max{|/i| < 1 : /i G a(/^(^, A))}, /i« = min{|/i| > 1 : /i G a{fxi^, A))}. 

Our second result shows that the constants cx, play an important role in 
the behavior of the bifurcation function gs{T, A). 

Theorem 5 Let the assumptions of Theorem IJ\ hold. Then the functions x%^s, 
gs have the following properties 

(i) 

Xw.,i3s{f3T, A) = f3xz^s{r, A), where I{f3s) = I{s) - 1, (21) 
g^s{f3r,X) = /3gs{T,X). (22) 
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(ii) For some constants C > 0, a > 0, independent of s,N and i G I{s) 

|(7,(r,A),-(cA(A-A) + c,r|)| <C((A-A)2 + |A-A|||r|U + ||r||L + e-"^/'). 

(23) 

Remark 6 In order to apply unfolding theory to the bifurcation equations ( fT6l) 
one also needs estimates of derivatives in (123|) . Though the functions (?s(-, ■) 
come out smoothly from Theorem [73, estimating their derivatives seems to be 
quite involved (cf. the proof of (l23l) in Section\M) and has not yet been done. 

If we consider a homoclinic symbol s E VLn with K = card(/(s)) < oo humps, 
then the theorem shows that the bifurcation equations are small perturbations 
of a set of K identical turning point equations 

= CA(A-A) + c,r2, iel{s). 

If ca, 7^ one can shift A to zero and scale A and such that one obtains a 
set of hilltop bifurcations of order K, cf . [6l Ch. IX, §3] , 

= A-r/, £el{s). (24) 

For two-humped orbits the set I{s) contains two elements and the solution curves 
of fl24|) are shown in Figure [71 This case will be crucial for understanding the 
global behavior of homoclinic curves in the next sections. 



3 Homoclinic orbits and their continuation 

A typical example, which plays the role of a normal form for quadratic two- 
dimensional mappings, is the famous Henon map, cf. pTl [T9l [5], [10] which is 
defined as 

/(x,A)=(^+f_^^-^^?^ 

This map has fixed points 

= (l-MA)) ' ""^^^ = ^ ^ V1 + 25X 

and for A = 0.35 a transversal homoclinic orbit xs(A) w.r.t. the fixed point 
exists, satisfying Assumption A4. 

For numerical computations, we approximate an infinite homoclinic orbit 
X'z{X) by a finite orbit segment xj, where J = [n_,n^] fl Z. The segment is 
determined as a zero of the boundary value operator 



Tj{xj,X) 



Xn+i- fixnA), n = n_,...,n^ 

b{Xn_ , Xfi, ) 
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Here b : R^'^ — )■ defines a boundary condition, for example 



-'per 



X 



or 



^proj 



fx 








in case of periodic and projection boundary conditions, where Bg and yield 
linear approximations of the stable and the unstable manifold. Due to our hy- 
perbolicity assumption, rj(-,A) has for J sufficiently large a unique zero in a 
neighborhood of the exact solution. Moreover approximation errors decay at an 
exponential rate that depends on the type of boundary condition, cf. |4]. 

For Henon's map, we solve the corresponding boundary value problem, obtain 
in this way an approximation of xs(A) and continue this orbit w.r.t. the parameter 
A, using the method of pseudo arclength continuation, cf. [EllIllH]. In Figured! 

we plot the amplitude of these orbits amp(xj(A)) := {J2nej W^nW — 'C+(A)|P)^ 
versus the parameter. 




ro,i = AO'i) 



0.35 0.4 0.45 0.5 



Figure 1: Continuation of homoclinic Henon orbits. At the parameter A = 0.35 
four distinct orbits x\, % G {0, ... ,3} exist that turn into each other via left or 
right turning points. 

At the value A = 0.35 four distinct homoclinic orbits occur that we denote 
by 2 G {0, . . . , 3}. We choose their index by following the order given by the 
continuation routine. The orbit is shown in Figure [2] together with parts of 
the stable and the unstable manifold of the fixed point ^+(A). The enlargement 
in this figure shows where the four orbits lie in the intersection of manifolds. 

At each turning point, two orbits collide; with r and i, we distinguish right 
and left turning points. Figure E] illustrates intersections of stable and unstable 
manifolds at these four turning points. 

Errors of turning point calculations for finite approximations of homoclinic 
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Figure 2: Primary homoclinic orbit x\ w.r.t. the fixed point (^+, and parts of the 
stable manifold (green) and the unstable manifold (red). The enlargement shows 
the intersections of manifolds that lead to the four homoclinic orbits in FigureUl 

orbits decay exponentially fast w.r.t. the length of the computed orbit segment, 
cf. [ini Theorem 5.1.1]. 

4 Connected components of multi-humped or- 
bits 

For Henon's map, we find four distinct transversal homoclinic orbits x^, s G 
{0, . . . ,3} at A = 0.35 and we identify with its symbol s. Note that the 

orbits 0, 1, 2, 3, pass into each other via left (L) and right (R) turning points tq i, 
^1,2, ■'^2,3, -^3,0, see Figure m The graph in Figure H] gives an alternative illustration 
of these transitions. 



L 




Figure 4: Transition graph for one-humped orbits. 

For the construction of an n-humped orbit, we choose a sufficiently long in- 
terval J = [n_, n-|_] around zero and a sequence s G 5" := {0, . . . , 3}". We define 
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Figure 3: Intersections of stable and unstable manifolds at the four turning points 
in the cutout region from Figure [B 

the pseudo orbit 

see Figure [51 Since the collection of single orbits x^, r G {0, ...,3} forms a 
hyperbolic set, the Shadowing-Lemma, cf. [23l [21] shows that the pseudo orbit 
lies close to a true n-humped /-orbit which we denote by In iS" 

there are 4" different symbols and thus we expect to find 4" different n-humped 
orbits We identify these orbits with their symbol. 

Note that the construction of pseudo orbits in (125|) slightly differs from f[T2|) . 
where we add up shifted orbits. With both approaches, we expect to find the 
same shadowing orbit for sufficiently large intervals J. 
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Xn 



Xn 



a:s[123] 



•^a • • • 



• • • • 



Figure 5: Construction of multi-humped orbits. 

Given two symbols s, s G iS", we analyze whether the n-humped orbits 
and x%{s) can turn into each other via continuation. 
Let us first look at the two-humped case. 

4.1 Bifurcation of two-humped orbits 

The continuation of two-humped orbits exhibits three closed curves of homoclinic 
orbits, cf. Figure ini and at the parameter value A, 16 different homoclinic orbits 
s G iS^ exist. 

amp 



3.8r 3.8r 3.8 




S.1 0.2 0.3 0.4 §.1 0.2 0.3 0.4 §.1 0.2 0.3 0.4 



A 

Figure 6: Continuation of two-humped orbits of length n„ = —20, = 21. 

One observes that at each turning point in Figure [6|, exactly one component 
of the symbol changes. For example, the symbol (1, 1) changes at a left turning 
point into the symbol (2, 1). Due to the combination of humps with finite length. 
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a perturbed hilltop bifurcation decides, whether (1,1) bifurcates into the symbol 
(2, 1) or into (1,2), see Figure [71 

For two-humped orbits, the system ( l24l) is a set of two equations in three 
variables A,ri,r2, called the hilltop normal form, cf. [6] 



A 



1 1 



A 



(26) 



Figure [7] (left) shows the solution curves of (126|) while the red curves in Figure 
[7] (right) indicate the generic solution picture of a perturbed equation. Here 
we neglect more detailed bifurcation diagrams which take into account small 
hysteresis effects w.r.t. the parameter A, see [6l Ch. IX, §3] for the unfolding 
theory. 




Figure 7: Unperturbed (left) and perturbed hilltop-bifurcation (right) at the turn- 
ing point ii^2- 



4.2 Connected components and equivalent symbols 

Homoclinic orbits that lie on a common closed curve define a connected compo- 
nent of 

n := {(z/s. A) e i^{R') X R : = /(|/„, A) Vn G lim y„ = ^(A)}. 

More precisely, let s E and denote by C{s) C the connected component 
that satisfies (a:s(s),A) G C(s). 

Then, we obtain an equivalence relation by identifying two sequences s, s G 
S^, if the corresponding orbits lie in the same component i.e. 

s^s ^ (xs(s),A) G C(s). (27) 

In the following, we discuss how to find these equivalence classes. Particularly, 
we show under some generic assumptions that each equivalence class has at least 
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four elements, and for n-humped orbits it turns out that one class has at least 
An elements. 

For this task, we introduce a labeled graph G with vertices s G 5". Two 
vertices s and s E are connected with an L or i?-edge, if ^^(s) bifurcates into 
Xzis) via a left or right turning point. Since we do not know the effect of the 
perturbed hilltop bifurcation a priori, we put an edge, if the transition is possible 
for at least one perturbation. For example, the vertices (1, 1) and (2, 1) as well 
as (1,1) and (1,2) are connected with L-edges in case n = 2, see Section [4.11 
Precise rules for constructing this graph are stated in Section 14.31 

Our hypothesis is that the desired equivalence classes correspond to a special 
decomposition of this graph into disjoint LR-cjdes. 

In case of one-humped orbits, the only Li?-cycle is 01230, see Figure HI Con- 
sequently, all symbols lie in the same equivalence class, which matches the fact 
that all one-humped orbits lie on the same closed curve and thus, in the same 
connected component of Ji. 

4.3 Graph structure of homoclinic network 

In this section, we give precise rules for defining the labeled graph G which we 
identify with its adjacency tensor with entries R and L. 

First, we assume that only one of the n humps can turn into a neighboring 
hump at a turning point. 

Rl There is no edge from s G 5" to s G 5" if 

n 

s = s or ll-s - s\\i = y^^djsj, Sj) > 2, 

1=1 

where d is the distance on the cycle 01230. 

Now let s,s E 5" and assume \\s — s\\i = 1, then there exists a unique j such 
that Sj 7^ Sj. 

From A^'^ > A°'^ we conclude that the right transition at r2,3 can only occur 
if the orbit contains no and no 1 hump. Therefore, we define i?-edges in G 
according to the following rule. 

R2 G{s, s) = Rii 

• {s„s,} = {0,1}, 

• {-^j) ^j} — {2? 3} and Si G {2, 3} for alH = 1, . . . , n. 

Similarly from A^'^ > A^'° we conclude that the left transition at £3,0 can only 
occur if the orbit contains no 1 and no 2 hump. Our rules for L-edges are: 

R3 G{s, s) = L if 
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, n. 



• {sj, Sj} = {0, 3} and Sj G {0, 3} for all i = I, 

We expect a connected component to correspond to an LR-cjcle in this graph, 
i.e. a cycle on which L and i?-edges alternate. A precise statement of our hy- 
pothesis is as follows. 



Hypothesis 7 The connected components of n-humped orbits and thus the equiv- 
alence classes fl27|) are in one to one correspondence to a partition of G into 
disjoint LR-cycles. 

In case n = 2, the L and i?-edges are shown in the left and center picture of 
Figure [HI The right diagram additionally shows the Li?-cycles that correspond 
to the connected components from Figure [61 

Similar diagrams for n = 3 are shown in Figure [9] 



Figure 8: L-edges of G (left) and R-edges (center) for two-humped orbits. The 
cycles in the right figure correspond to the closed curves that are computed nu- 
merically in Figure\^ 




000 



Figure 9: L-edges of G (left) and R-edges (center) for three-humped orbits. The 
green lines in the right figure show the cycles that are computed numerically. 
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We continued n-humped orbits numerically for Henon's map up to n = 5. 
Table [1] summarizes the number of cycles and their lengths found in the compu- 
tation. 



n 


le 

4 


ngth of cycle 
8 12 16 


20 


1 


1 




2 


2 


1 




3 


9 


2 


1 




4 


45 


3 


3 


1 




5 


205 


6 


6 


4 


1 



Table 1: Continuation of n-humped orbits - length of occurring cycles in numer- 
ical experiments. 

These experiments for n-humped orbits of Henon's map suggest that the 
length of occurring cycles is a multiple of 4. Furthermore one cycle exists of 
at least length 4n. 

Theorem 8 Fix G N and assume that Hypothesis \7\ holds true. Then, all 
n-humped orbits lie on cycles whose length is a multiple of 4. 

For the specific symbols s° = (0, . . . , 0) and = (2, . . . , 2) G 5", the corre- 
sponding orbits and x^is"^) He on a common cycle of at least length An. 

Remark 9 Table U\ shows that there is exactly one orbit of length 4n and all 
other orbits are shorter. Hence, the orbit of length An contains s'^ and . 

Proof: By assuming Hypothesis [7] we see that it suffices to analyze Li?-cycles of 
the graph G. More precisely, we prove Theorem [8] along the following steps. 

(i) Each LR-cjcle in G has length 4m, m > 1, m G N. 

(ii) There exists an LR-cjc\e from s° to of length An and each Li?-cycle that 
contains s° and has at least length An. 

(iii) Each LR-cjde that contains also contains s^. 

(i) Let the vertices , . . . , v'^ , v"^~^^ = form an Li?-cycle in G. Fix j G 
{1, . . . ,n} and let Ij = #{z g {1, . . . , m} : ^ vj+S G{v\v'+^) = L} be 
the number of L-edges for which the corresponding vertices only differ in 
the jth component. 

From Rl it follows that Ij is an even number and maXj=i_,,,^„{Zj} > 2, 
otherwise the cycle cannot be closed in the jth component. Furthermore, 
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an Li?-cycle has the same number of L and i?-edges. Thus, the cycle has 
length 

n n n 

m = = 4^7^ = 4p with P = ^7:>^- 

i=l i=l i=l 

(ii) We explicitly construct an Li?-cycle in G from s*^ to s^: (0 ■ ■ ■ 0) i— (10 ■ ■ ■ 0) 
K> (20 ■■■0) ^ (210 ■■■0) ^ (220 ■■■0) ■ ■ ■ (2 ■ ■ ■ 2) ^ (32 ■ ■ ■ 2) ^ 
(312 ■■■2) h-> (302 ■■■2) h-> (3012 ■■■2) ^ (3002 ■■■2) ^ (30- --O) h> 
(0 ■ ■ ■ 0) which has length An. Note that the distance from s° to on 
the full quadratic grid is 2n and consequently, each cycle containing these 
two points has at least length An. 

(iii) For proving that each Li?-cycle with s° also contains s^, assume that an 
LR-cjcle exists that contains but not This cycle lies in the subgraph 
that we obtain by deleting the vertex s° with its corresponding edges. 

Note that the Li?-cycles start with an L and end with an i?-edge, and it 
follows from R2 that all i?-edges that start at end in 

Vs := {s e S'' : 3j e{l,...,n}: Sj = 3}. 

We define the graph G by deleting the i?-edges of from the remaining 
graph. Figure [TO] illustrates this construction in case n = 2. 
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Figure 10: Transition graph in case n = 2. L-edges and R-edges are plotted in 
red and black, respectively. Vertices and edges that are deleted in the proof of 
Theorem\E are marked by green crosses. 

If G breaks into two components V3 and G\V3, then we get a contradiction 
to the above assumption and an LR-cjcle in the original graph G that 
contains but not cannot exist. 



17 



To finish the proof, we show that an Li?-path in G from to V3 does not 
exist. 

From we cannot go directly via an i?-edge to V3, since the corresponding 
edges are deleted in G. Thus without loss of generality, we get: (2 ■ ■ ■ 2) i— )■ 
(12 ■ ■ - 2) I—)- (02 ■ ■ - 2). Denote by t;™ the mth vertex on this path. A 3- 
component can only be achieved by the left transition £3 or by the right 
transition r2,3, see R2 and R3. 

If a j exists with t;™ G {0, 1}, then the r2,3 transition is impossible by R2. 

If a j exists with v™" G {1,2}, then the transition is impossible by R3. 

Thus, we only obtain a 3 component via the vertex s° = (0 . . . 0) which is 
deleted in G. 



5 Bifurcation analysis near homoclinic tangen- 
cies 

In this section we prove the main Theorem [1] and Theorem [5t^i) by using an 
existence and uniqueness result for a suitable operator equation in spaces of 
bounded sequences. In the following we use the notation Bp{x) and Bp = Bp{0) 
to denote closed balls of radius p in some Banach space. 



5.1 The operator equation 

First recall the operator F : ^°°(R^) x R — )■ £°°(R'^) from fl2U]) and the normal- 
ization A = and .^(A) = for A close to A, see A2, A3 in Section [2J Then for 
any s G fi^v define the operator 

Gs : £~(R'=) X £°°(s) X i°°{s) X R ^ i^{R'') x £°°(s) 

by 

Here p^, are defined in ( |T2l) . ( |T4l) and w{s, g) is given by (recall W'z from (fTTll ) 

w{s,g) = J2 9eP~^w^, 9 e £°°(s). 

Our aim is to derive the functions gs in Theorem [U by solving 

Gsix^,g,T,X) = (29) 

for ||t||oo, |A| sufficiently small and for all s G Qn- More precisely, we prove in 
Section [6] the following Reduction Theorem. 
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Theorem 10 There exist constants Cq, px, Pg, Pr, Px > and a number Nq & ¥i 
such that for all N > Nq and for all s & Qn the following statements hold. For 
all T G Bp^, X G Bp^ the system ( !29l) has a unique solution 

g = gsir, A) G Bp^ C £^{s), = Xs,,(r, A) G Bp^ C £°°(R'). 

Moreover, the following estimate is satisfied: 

max(||5( - ^lloo, \\xz -Xb||oo) < Co\\Gsixz, g,T, \) - Gs{x^,g,T, \)\\ (30) 

for all g,g e Bp^, G Bp^, r G Bp^, \ E Bp^. 

5.2 Preparatory Lemmata 

In order to construct the neighborhoods U and A in Theorem [1] we need several 
lemmata. 

Lemma 11 There exists Ni G N, such that for all N > Ni, s E the linear 
system 

J2 (/3"%,/3"%)rfc = r,, £ G /(s), r G (31) 
fce/{s) 

has a unique solution r G £°°(s) anc? 

Iklloo < 2||r|U. (32) 
Proof: Rewrite (13T|) as fixed point equation 

r = Pr + r, (Pr), = - ^ {^'u^, T^u^H 



and note 



|^r||oo < Iklloo sup ^ K^iS,^^ ''U'z) 



Ce 



^ 1 - e-"^ 



-ciiV- 



Thus, we choose A''! such that -^_^_^jv^ < |. Then P is contractive and (132|) 
follows. ■ 

Lemma 12 There exist N2 E ¥i , C2 > such that 

WFip^is) + v^is, r), A)|U < C2(|A| + ||r||L + e'^^/^) (33) 

/or a// iV > iVs, s G n^, |A| < 1, r G Si C £°°(s). 
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Proof: We estimate 

Pn+l{s) + Vn+l{s, r) - f{Pn{s) + Vn{s, t), A) 
= ^ {Xn+l-l + nUn+l-l) - f{Pn{s) + t) , 0) + C( | A | ) 

iei{s) 

-/( J] ^ r,u,_,,0) +0(|A|) 

^e/(s) £e/(s) 

-/( X a;„_,,0)-/4 X ^n-j,0) J2 nun^i + 0{\X\ + \\r\ 



' 1 



(34) 



For n E li, choose i G I{s) such that |r2 — ^| < \n — i\ for all £ G /(s). Then 
|n - ^1 > f for alH 7^ £ and, therefore, by ©, 

II X /(^n-^0)-/( X ^n-^,0) 

< II X /(a;n-£,0)|| + ||/(x„_,-,0)-/( X ^n-£,0) 

< II X a;„+i„,||+L|| X 3;„„,|| < Cse-"^/^ 

££i{s),ej^e ££i{s),e^e 

In a similar way, 

^G/(s) j£l(s) £elis) 

Combining these estimates, we obtain (l33l) . 
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Lemma 13 Assume Al, A2 and let x% he a homoclinic /(■, 0)-orbit with respect 
to the hyperbolic fixed point 0. Then, there exist zero neighborhoods C R*^, 
A3 C Aq and constants N-^,nQ,a > 0, C3 > 1 such that the following statement 
holds for all K > N^, —n_,n+ > Hq: 

If Xn+i = f{xn,X) for n e J := [0,K - 1], A G A3, and if Xn G U3 for all 
n E J := [0,K] then we have the estimate 

sup \\xj-Xn_+i-K+j-Xn+-i+j\\ < C3 (\\xk " ^n_+i|| + ||a;o " + |A| + e""^ ) . 

i6J ^ ^ 



Furthermore we obtain in case K = oo: 



sup \\Xj - Xn+-l+j\\ < CsiWxo - Xn+-l\\ + |A|). 
j>0 



(36) 
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Proof: Consider the pseudo-orbit pj := x„_^_i+j + Xn_+i-K+j, j ^ J which is a 
zero of the boundary value operator 

V bK{yo,yK) J 

at A = 0, where 

Pn ■ = 0) + /(a^n.+i-x+n, 0) — 0), n&J, 

Here Pg and P„ are the stable and unstable projectors of the fixed point 0. 
Let b be the bound from Theorem [2T] for the difference equation 

Un+l = (/x(0, 0) + Bn)Un, = fx{Pn, A) - /^.(O, 0), UEJ. 

For sufficiently large — n., ?7,+ > no and A G A3 sufficiently small, we get ||-Bn|| < h 
for all n E J . Consequently 

Un+l = fn{Pn, >^)Un, U G J 

has an exponential dichotomy on J with projectors P^,P^ and an exponential 
rate a that is independent of n_, n+, A and K. 

As in the proof of [12], Theorem 4], we show that for A G A3, n_, n+ > Uq and 
K > N3 we have a uniform bound 

\\DiTjipj,\)-'\\^<a-\ (37) 

In order to see this, consider the inhomogeneous difference equation 

Un+l - fxiPnA)Un = r^, U = 0, . . . K - 1, (38) 

PsUq + PuUk = 7- (39) 

Denote by $ the solution operator of the homogeneous equation and let G be the 
corresponding Green's function, cf. ( |87l) . The general solution of ( |38l) is given by 

Un = ^{n,0)v + ^G{n,m + l)rm, (40) 

where 

= + $(0, K)v+, G 7^(Po'), V+ G 7^(P^). 

Inserting ( HOl) into ( l39l) . it remains to solve 

P,(t;_ + <D(0, ir)t;+) + P„(<l>(i^, 0)t;„ + = P, 
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with 

This finite-dimensional system has a unique solution for K > sufficiently large 
since ||Ps -Pg || and ||P„ -P^|| -)■ as i^' -> oo. Therefore, the system fl55]l . 
fl5^ also has a unique solution mj for K large and the dichotomy estimates lead 
to a bound 

i.e. da?]) holds. 

We apply Theorem [T^ with the space Y = Cf = {(2/n)n6J • Un ^ K-'^} of finite 
sequences and with Z = x. R*', both endowed with the sup- norm. We take 
Vo = Pj and use uniform data for all A G A3. For 6 sufficiently small we have 

WDiTjixj, A) - DiTji^pj, A)||oo < ^ for all sequences xj e Bs{pj), A G A3, 

and by choosing the neighborhood A3 sufficiently small we get 

||rj(]?J, A)||oo = SUp\\pn+l- fiPnA) - Pn\\ + \\bKiPO,PK)\\ 

neJ 

= sup \\pn+l - f{Pn, A) + /(p„, 0) 
nGJ 

0) - f{Xn_+l~K+n, 0)11 

= snp \\ f {pn, 0) - f{pn,X)\\ < ^6 for A e A3. 

Theorem [19] applies to A G A3 with uniform data, and it follows from f lllQp with 
some constant C3 > that 

\\xj-yj\\oo < C3||rj(a;j, A)-rj(?/j, A)||oo for all xj,yj G Bs{pj),\ G A3. (41) 

From ([7]) we find a number such that x„ G 5^(0) for all |?t,| > uq and also 
Pn £ P*(0),n G J for all — n_,n+ > uq. Then we take ?73 := Bs_{0) as our 
neighborhood and note that ( l4Tll holds for any two sequences xj,yj in f/3. For 
n & J and A G A3 it follows that 

/bn, A) - Pn+1 

= f{^n+-l+n + 2;n_+l-_X+n; A) — — Xn_+l-K+n+l 

/'('^n+ — 1+71 ~l~ "^n- +l_X+n, 0) - f{Xn+-l+n, 0) - f{Xn_+l-K+n, 0) + C(|A|) 
-fiXn_+l-K+n, 0) 0{\\Xn_+l^K+n\\) for < Tl < f 

0) + C(||x„+_i+„||) ior f<n<K 

C(e-"f + |A|). 
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Now let xj be a sequence in such that Xn+i = f{xn, A) for all n & J, and 
some A G A3. Then 

\\xj-pj\\oo < c||rj(xj, A) - rj(pj, A)||oo 

< C(^SUp ||/(p„,, A) -Pn+l\\ + \\bK{xo,XK) 



n<=J 



< C e-"- + |A| + 



Ps{xo-po) 
Pu{xk -Pk, 

< C (e~°^ + ||a;o - Poll + W^k - Pk\\ + |A| 



< C3 (^e""^ + ||a;o - ^n+-i|| + I^a' - ^n_+i|| + |A| 
In case K = 00, one uses the operator 



and it turns out that Dif ]n(p]n, A) with pj = Xn+-i+j has a uniformly bounded 
inverse for A G A3 and sufficiently large — n_, n+. Then the estimate (l36l) follows 
immediately. 



5.3 Proof of Main Theorem 

Let us first prove assertion (i) in Theorem [H 
Step 1: (Construction of neighborhoods U,A) 
In the following Ai D A2 D . . . will denote shrinking neighborhoods of 0. Let 
Pg, Px > he given by Theorem [TO] and note that we can decrease Pt,P\ with- 
out changing the assertion of Theorem [TOl Introduce the constants (cf. ([7]) and 
Lemma [111 [ED 

a, = e-^ C, = Cs + -^, C, = Ce^^, C, = CoC2 + -^. (42) 

1 — a* 1 — a* 1 — a* 

Let Pt- > be such that 

C,Pr < f . (43) 

By Lemma [T^ we can choose a ball B^/.^ C f/3 and numbers n+, — n_ > uq such 
that Xn G -Beg for all n > n_|_ — 1 and n < n_ + 1. It is well known that the only 
full orbit in a small neighborhood of a hyperbolic fixed point is the fixed point 
itself. That is, we can assume w.l.o.g. that f/3, A3 satisfy 

Vn+i = fiVn, A), Vn & U^{n G Z), A G A3 ^ = for all n G (44) 
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The set /C = {0} U {xn : n < n_ or n > n+} is compact and satisfies 

/(/C,0) C/CU{x„_+i}. 
Thus we find a.n e < Eq and A4 C A3 such that the following properties hold 



Vo := 5, U U C B2e, C U- 



3; 



the balls 



Bj := B^{xn^+j), j = 1, . . . , X := n+ - 1 - n_ 
are mutually disjoint, 

/(yo,A4)n5, = for j = 2,...,x, 



2^5(2^3 + 1)£ < ^ and {2Cs + l)e<^ 



see Figure [TT] 



Bi 




Bo 





B-^-i 




n+ 





(45) 

(46) 

(47) 
(48) 



•^j-j I — 1 



Figure 11: Construction of neighborhoods. 
Next we take A^4 > max(A^i, A^2, ^3) (see Lemmata [TTl [T2| [T3|) such that 
2C5C4e-"^^ < ^, C4e-"^^<^, {2{C, + Cs) + 18CsC,)e-'^'^- < (49) 
We can also find A5 C A4 such that for all A G A5 

2C5C3|A|<^, C3|A|<^, \M<Px, 6C3|A|(3CoC2 + l)<£. (50) 



Finally, we define 



N = X + N^, where iV* = 2A^4. 



(51) 
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Then we choose Ag C A5 such that the following settings define neighborhoods 
Vj of Xn_+j,j = X, . . . , 1 recursively (cf. Figure [TTj) : 



n=l ^ ' 

V,=B,nr\V,+i,Ae), for j = x- 

Here we use the notation /""(V^-, Ag) = {x : /"(x, A) G Vj for all A G Ag}. 

With these settings we consider the maximal invariant set M(f/, A), cf. (j3]), 
that belongs to 

U=[jVj, A = Ag. (53) 

j=0 

Let us note that our construction (146|) . (147|) . ( |52|) . ( 153|) implies the following three 
assertions for any /(-, A)-orbit ?/a C t/, A G A 

l/n e Vo, ^Vq^ yn+1 G Vi, (54) 

Vn G V,- for some l<j<x-l^ G (55) 
Vn e ^ Vn+i e Vo for £ = l,...,iV, + l. (56) 
Step 2: (Construction of symbolic sequence s) 

For some A G A consider an orbit of ([I]) that lies in U. If it lies in Vq then 
Un = 0,n E hj (jH]) and we set s = G I^at. Otherwise we have yn ^ Vq for 
some fi G We show that 

Hy^) ■={neZ: yn G V^} 

is nonempty and that there is a unique s G Vt^ such that /(?/s) = -^(-s), see ([9]). 
By our assumptions we have yn G Vj^ for some jo ^ {1; • • • ? If Jo = 1 then 
-^(z/z) 7^ whereas in case jo > 2 we obtain yn-m £ Vjo-m, ""^ = 0, . . . , jo — 1 by 
induction from ( 155|) and the fact that the Vj are mutually disjoint. Therefore, we 
have £ := n — jo + 1 G I^yw.)- Moreover, from (l54l) and (l56l) we obtain 

|/^-+^. G i = 0, ... X - 1, y^-,^^. G Vo, j = X, . . . , X + iV,. (57) 

This shows that the difference of two consecutive indices i < i in lilt's) is at least 
= X + A^^,. Therefore /(|/s) belongs to Z(A^) (cf. (fTOl) ) and there is a unique 

sequence s G I^at such that /(ys) = -^(-s)- 

The relations fl57j) hold whenever £ G /(s). By P6l) and (152!) this gives us the 

estimates 

\\yp+t+u - XpW < e for £ G /(s),p = n_ + 1, . . . ,n+ - 1 (58) 
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for the index v = — n_ — 1 (cf. ( IT5|) ). From this we will derive the estimate 
Tn := WVn+u - Yl ^ ^^46-""^^ + C3IAI + (2C3 + l)e for neZ. (59) 

mel{s) 

Consider first indices n = p + £ with i G I{s) and p = + 1, . . . — 1. Then 
with ( 15T1) we find 



(60) 



1 - «* 

Next consider two consecutive indices i < i in I{s) and n = p + £ for p = 
— l,...,p = ^ — £ + + For these indices we get 

{V^, for p = n+ - 1, 
Vi, for p = ^, 
Vo, otherwise, 

and we can apply Lemma fT3l to this sequence in place of xq, . . . , xk, where K = 
i-i + 2-{n+-n^)=£-i+l-H>N, + l = 2Ni + l. With dSHD this yields 
the estimate 

sup \\ye-^u+n+-l+j ~ Xn++j-l ~ Xn_+j-K+l\\ 
j=0,...,K 

< Cs {\\yi+^+n+-i - Xn+-i\\ + \\yi+u+n+-i+K " Sn_+i|| + |A| + a''^^'^) (61) 
<C3(25+|A|+«f^). 

Finally, we use this to estimate for ra = p + £ and p = — 1, . . . ,p 

Tn ^ ll^/n+i/ ^n—£ + ^ ^ H^^n— m||- 

in&I{s),m=/=e,i 

The first term is handled by ( 16T1) . Further note that 

E ll^n-.n|| <CeJ]af*<f^. (62) 

- ^ ^ I - a* 

e<mei{s) J>i 
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and the same estimate holds for £ > m G /(s). Collecting estimates (160|) to (162|) 
we arrive at (159|) : 



2r' 

T„ < (2C3 + 1)6 + C3IAI + Caaf + < C4e"-^^ + C3IAI + (2C3 + l)e. 

1 — a* 

Finally, we note that (159|) also holds in case £ is the smallest index in J(s) or i is 
the largest index in /(s), respectively. Then one repeats the previous arguments 
with the formal setting i = —00 resp. i = 00 and uses the corresponding one- 
sided version of Lemma [121 

Step 3: (Construction and estimate of r and a;^) 

We want to find G i°°{R''),T E i°^{s) such that ([29]) holds with g = 0. The 
second term of the operator G in f[28|) vanishes provided we solve 



J] (/3-%,/3-'"Ms)r„ = -Ps(s)), ^ G /(s), (63) 

and set 

a^z = /^''yz -Pz(s) - fz(s,r). (64) 

By Lemma [TT] the linear system ([631) has a unique solution r G ^°°(s) which 
satisfies bv (|59|), ([7[), ([12|) 



klloo < 2||ms||£1 sup - ^ 



mGJ (s) 

< 2C}-^ (C4e-"^^ + C3IAI + (2C3 + l)e) . 
1 - a* 

Using ([IH}, ([IHD, ([50D and CeYzf; < C's we end up with ||r||oo < Pr- Next we 
estimate from ([64ll . Using ([7]) it is easy to show that 

II T r%||oo<a^-^ = C5. (65) 
Therefore, when using (15 9 1) again, we find 

Ikzlloo < Wf^^yw, -pz(s)||oo + C^Pr 

< C4e-"^* + C3IAI + (2C3 + l)e + C^Pr. 

The estimates from fll3l) . psj) . fll9l) . f[50|) then guarantee ||xa||oo < Ps- 

Therefore, we know that Gs(xa,0,r, A) = and the tuple (xs,0,r. A) lies 
in the balls in which (129 p has a unique solution. By uniqueness we conclude 
gs{T, A) = and a;^ = xs,s(t. A) for A G A. Moreover, by the defining equations 
( 163|) and ( 164]) we obtain that equality (fTSjl holds. 
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step 4: (Proof of Theorem [T] (ii)) 
The radius r^- will be taken such that 

18C3(CoC2r2 + C,r,) < e. (66) 

Let r G B^^ C ^°°(s), A G A satisfy Qsij, A) = for some s eVLn and let be 
given as in Theorem [101 Then clearly, the sequence 

Vtl = a;z,,(r, A) + ps(s) + v^^s, r) 

is an orbit of ([1]). It remains to show that yn E U for all n G Application of 
(j30|) in Theorem [10] and of Lemma [12] yields the estimate 



||a^z,s(^, A)||oo < Co\\Gs{x'z,s{r, X),gs{r, A),r, A) - G,(0,0,r, A)||oo 

= Co||F(ps(s)+t;s(s,r),A)|U (67) 
<CoC2(|A| + ||r||L + e-°^/2). 

From f[65]) we find 

\Ms,T)\\oo<C5rr. (68) 
We estimate the distance of to the centers of the balls Bj in U by showing for 

U = —71- — 1 

WVp+i+u - XpW < for p = n_ + 1, . . . ,n+ - 1,£ G /(s). (69) 

6C3 

Note that the right-hand side is less equal e since we chose C3 > 1 in Lemma [T3l 
Using ([67]), ([68D and ([60]) we obtain for I G I{s) 

\\yp+i+y-Xp\\ = \\xp+i^s{T,\)+ ^ + fp+£(s,r)|| 

< CoC,{\X\ + llrllL + e-"^/^) + Csr. + -^af^ ^^^^ 

1 - a* 

< C0C2IAI + (CoCsr. + C5)r. + C,a^'. 



Conditions ( \i9\i . ( 150]) and ( 166]) guarantee that ( 169]) is satisfied. 

Next we consider two consecutive indices £ < i in I{s) and n = p + i for 
p = n+ — 1, . . . ,p = i — i + n_ + 1. In the first step we show that yp+£+u ^ U3. 
Using p > — 1 and p + i — i<n_ + l and fl62]) we estimate similar to f lTO]) 



A)|| + || + Yl \\^P+^~rn\\ + \\Vp+i{s,T) 

or 

< C,Cm + llrllL + e-"^/2) + Csr. + -^af^ + 25o 

1 - a* 



< C0C2IAI + (CoCsr. + C5)r. + Cgaf* + 25, 
6C,3 



28 







Since i^s^Q C U3 this proves our assertion. Now we can invoke Lemma [13] and 
find as in (pTl) 

sup Wvi+^+n+^l+j - Xn++j-l - Xn_+j-K+l\\ < C3 ( + |A| + Ctf M . (71) 

j=0,...,K V'3<-'3 / 

For j = 0, . . . , A^4 we have + j — K + 1 < —N4 + n_ and hence 

II ■ ,\\ < C pa{n-+N4-K+l) ^ ^ N4, 



while for j = + 1, . . . K we have + j — 1 > ri.!. + A'4 and hence 

IK?i++j-l|| ^ '-^e«* ^ • 

Combining this with ( 17T|1 we find 



- Xpll < I + C3IAI + (C3 + Ce)af^ p = n+ - 1, . . . ,n+ - 1 + Ar4, 
- Xp_l_^_f,\\ < I + C'alAI + (C3 + Ce)afS p = n+ + A/'4, . . . ,p. 

We have arranged the constants in ( l49l) and ( !50l) such that the right hand side is 
bounded by e. 

For the final step we note that we just have shown (cf. ( H5l) ) 



Kr, forp = n+-l, 
yp+e+^ e { Vl, for p = p = £ - £ + n_ + 1, 
Vq, for p = n+, ...,£ — £ + n„ 



for two consecutive indices i < i in I{s). On the other hand we know from 
that Up+t+y E Bs{xp),p = n_ + 1, . . . , n+ — 1. Since is an /(-, A) orbit we 
conclude by induction from the definition ( 152|) 

Vp+e+y e Vp_n- , for p = n+ - 1, . . . , n_ + 1. 

Therefore the sequence ?/s hes in U which proves our assertion. 

Step 5: (Proof of Theorem [5] (i)) The proof of ( 12T|) . ( l22i) is easily accomplished 
by noting the equivariance relations 

/(/3s) = I{s) - 1, pw.{(3s) = (3p'e{s), Vz{(3s, (3t) = (3v^{s, r), 
w{/3s,/3g) = I3w{s,g), F{f3xz, X) = f3F{xz, \), and 
= (/3'+V,/3xs) for iel{l3s) = l{s)-l. 

The assertion then follows by uniqueness from Theorem since neighborhoods 
are shift invariant as well. 

The proof of Theorem EJ^ii) will be deferred to the next section. 
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6 Proof of Reduction Theorem 



6.1 Nonlinear estimates 

According to ( !28l) the Frechet derivative of Gs w.r.t. X'^^.g is given by 



D(^^^g)Gs{xz,g,T,X){y^,h) 



'D^F{p'z{s) + + r), A)i/s + w{s, hf 

D^F{X'Z, \)y7L ={yn+l - fx{Xn, \)yn)n(i7L- 

(72) 

The key step in the proof of Theorem [10] will be a uniform bound for the inverse 
of D(,,3)G,(0, 0,0,0). 

Lemma 14 There exist constants C,N > such that for all N > N , s ^ Qn 
the operator D(^^ g)Gs{0, 0,0,0) is invertible and satisfies 



\y^\\oo + \\h\\oo < C'p(.,3)G,(0,0,0,0)(i/s,/i)||, E £°^(R'=),/i G 



s . 



Before proving Lemma [H] we finish the proof of Theorem [101 
Proof: Let and Lx be Lipschitz constants of the Jacobian fx{x. A) with respect 
to X and A in a compact ball that contains the homoclinic orbit in its interior. 
Then formula ( !72|) and the bound (165]) directly lead to the Lipschitz estimate 

|lA^,3)G's(x^,^i,ri, Ai) - D^x^g)Gs{xl,g2,T2,X2)\\ ^^^^ 
< L^iWxl^ - xlWoo + Csllri - Talloo) + ^a|Ai - A2I 

for all x\,x\ G Bp^, g\,g2 G , ti,T2 G Bp^, and Ai,A2 G Aq with Px,Pt taken 
sufficiently small. From Lemma [TJ] and Lemma [TS] we obtain that the operators 
D{x,g)Gs{0, 0, T, A) are invertible for r G Bp^, A G Bp^ provided we choose 

N>N and L^G^Pr + ^^aPa < i- 

Then we have 

\\D^x,g)G,iO, 0, T, Xy'W <2C, TEBp^,\E Bp^. 

Now we apply Theorem [T^ to every operator F = Gs{-, ■,t, X) in the spaces 
Y = Z = £^{R^) X £°°(s). Setting a = ^,yo = and taking L^px < ^ we find 
that condition f lll6p is satisfied with = ^ and S = px- Finally, we obtain from 
Lemma [H for all > A^2, 

||a(0,0,r,A)|| = \\F{p^{s)+v^{s,T),\)\\^ 
<C^2(|A| + ||r||L + e""^/^) 
<C2(pA + P^ + e""^/^). 
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Now we select Nq > max(iV, N2) and px, pr such that C2{px + Pr + e""*^") < ^. 
Then we find ||G,(0, 0, r, A)|| < (a - k)S = ^ for all > No, i.e. condition 
(I117p is satisfied. An application of Theorem |19] finishes the proof. ■ 

Remark 15 // instead of Theorem [7PI we use a Lipschitz inverse mapping the- 
orem with smooth parameters (of. 123. Appendix]) then it is easily seen that the 
solutions QsjX'z.s are smooth functions of t and A. 

Proof: (Theorem [5] (ii)) 

With ca, from (USD define h = h{T, A) G i°°{s) by 

hi = cxX + c^Tf, iel{s). 

The idea is to construct elements ?/s = A) G i°°(JR!') and 7 = 7(r, A) G 
such that the residual 

Gs{pw.{s) + 2/E + r), /i + 7, r, A) 

and 7 are of higher order than 0{\X\ + ||t||^). Then the assertion follows from 
fl30|) by comparing them to Xa,s(T, A), (^^(r, A). We find y^z, 7 by Taylor expansion 
of F (we abbreviate F° = -^^(Pzls), 0) etc.) 

+y^ + v^is, r), A) = F° + + v^{s, r)) + F°A 

+ ^F^Ay^ + v^is, T)f + Fl^{y^ + r))A + ^F^.A^ 

+ C((|A| + ||r|U + llyslloo)^). 

From the estimates (!34|) . (l35l) in the proof of Lemma [T2] we have 

||F°|U = 0(e-"^/^), ||F>^(s,r)|| = 0(e-"^/lr|U). 

In a similar way, using Lipschitz constants for fx and we find 

(Fl(t;s(s,r))2)„, = - ^'/xx(a:„-^,O)K_,)^ + 0(e-"^/2||r|U). 

Therefore, Taylor expansion of Gg yields 

GsipM + y7L + V7l{s, r), /i + 7, r. A) = D^^^g)G'l{yw„ 7) + (v^z, 0) 
+ 0(e-'^^/' + ||r|U|A| + \\r\\M\oo + A^ + WvAlo + IkllL), 

where 
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This suggests to define {y'z,^) by 

^(x,,)G'°(2/s,7) = -(<^s,0). (75) 

From this equation and Lemma [14] we have the estimate 

lbdloo + ||7l|oo<C'(|A| + ||r||^). (76) 

Taking the inner product of the first coordinate in ( 1751) with /S'^w^ji G I{s) and 
using flT9|l . leads to the improved estimate 

hiloo < Ce-^^/^dAI + llrllL). (77) 

By fl76l) . fl75|) the Taylor expansion fl74l) of assumes the form 

||a(Pa(s) + Z/s + ^s(s,r),/i + 7,r,A)|| = 0(e-°^/2^||r|U|A| + A2+||r||L). (78) 

With fl5U]) and fl77|) this leads us to the final result 

lkz,.(r,A)-ys(r, A)||oo + ||^7.(r, \)-h{r, A)||oo = 0(e-"^/' + ||r||oo|A|+A2+||r||L). 

■ 

6.2 Linear estimates 

In this subsection we prove Lemma O For any two integers ni < Ur and for any 
number a > consider the weight function 

Un = ujn{a,ni,nr) = min (e-^^"""'), i,e"("'-'^)) , n G (79) 

Note that w^; has a constant plateau of arbitrary width with exponentially de- 
caying tails on both sides. We also allow ni = — oo and = oo (but neither 
Til = rir = —oo nor ri/ = = oo), in which case has only one-sided decay or 
degenerates to the maximum norm if n; = — oo and 77,^ = 00. In the following 
we will suppress the dependence of the norm on the data ni,nr,a, but all our 
estimates will be uniform with respect to 

— 00 < rii < rir < 00 < a < ao < a, (80) 

where < ao < « is fixed. The following lemma shows that exponentially 
decaying kernels preserve the weight. 

Lemma 16 There exists a constant Kq, depending only on Oq, a, such that 
and for all weight functions satisfying (IHOj) . 
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Proof: We consider 

On — UJn 6 

■mew. 

only for n < ni and leave cases ni + 1 < n < rir and + 1 < n to the reader. 

ni-l 

r \ ^ a{ni—n)—a{n—l—m)+a{m—ni) , \ ^ a{ni—n)—a{m,+l—n)+a{m,—ni) 

On - ^ + / , 

m<r!.— 1 m=n 

rir — 1 

I ^ ^ ga(n; — n)— a(m+l— n) ^ ^ ^0(71; — n)— a(m+l— n)+a(nr — m) 
m=ni m>nr 

ni-1 

—Q-a Sr~^ g-(a+a){n-l-m) _|_ ^-a g-{a-a)(m-n) 



m<n—l m=n 
rir — l 

-j_g~("~")("i~")~" ^ ^ g-a(m-ni) _|_ g-(a-a){ni-n)+a(n;-nr-l) ^ g-{o+a)(m-nr) 
m=ni m>nr 



1 - e-(°+'^) V 1 - e-^"-"^) 1 - e- 

A suitable constant for all cases is Kq = 2 + e"° + -^_g-°°-ao) • ' 
With the weights from above we consider the Banach space 

C = {y^ e {R'f : \\y^\l = sup < oo}. 

Taking the exponent a from ([7]), f|T8l) we have the following result for the varia- 
tional equation IQ. 

Lemma 17 Let conditions Al - A3 and B4 hold. Then the linear system 

Vn+l - fxiXn, 0)yn + kWn = r„, 77, G Z 

{uw.,yz) = 7 

has a unique solution {y'z,h) G x R for every (rs,7) G x R. Moreover 
there is a constant C* such that for all weights (!79|) wi/i (IHOl) . 

<C*(||rsL + |7lhsL). (82) 
If ni<{] < Ur, then ( l82l) simplifies to 

+ |/i|<C*(||rsL + |7l). 

Proof: Let us abbreviate A„ = /^(x„,0), = (z/n+i - ^n2/n)nez and denote 
by $ the solution operator f lT^ of ffT^ . From [221 Section 2] (see also [15]) 
we obtain the following facts. Equation fll20p has an exponential dichotomy for 
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n > with data {K, a, P+^ P+") and forra < with data {K, a, P-^ P~") (see 
Definition [20]). Due to B4 we have 7^(Po+") n 7^(Po"'') = span{Mo} =: Yq and 
there exist decompositions 

n{P+') = Yo(BY+, 7^(Po-") = Fo©>^-, r+nr_ = {0}, (83) 

where 

fc, = rank(P+') = dim(F+) + l, (n > 0), 
fc„ = rank(p-") =dim(F_) + l, (n < 0), 
fc, = A; - A;,, P+^ + P+" = / (n > 0), P,"^ + P^ = / (n < 0). 

The operator C : ^""(R'^) £~(R'') is Fredholm of index with 

AfiC) = span{M4, 7^(£) = {r^ G ^°°(R'-) : {w^, r^) = 0}. (84) 

One can also choose the ranges of Pq^" and Pq~* such that, in addition to ( 183|) . 

R^ = © © F„ © , 7^(Po+") = r_ © Fi , 7^(Po"^) = F+ © n . (85) 

Here dimYi = 1 and one can take Yi = span{u7_i}. Following [221 Lemma 2.7] 
the general bounded solution of yn+i — AnUn = r.«, n > is given by 

y+ = <l>(n,0)77+ + 5^G+(n,m + l)r^, n > 0, r/+ G 7^(Po+^) (86) 

m>0 

with the Green's function defined as follows 

r(r,rr,\-l ^{n,m)P+\ < m < n, 
^+^'''"'^ ~ \ -<l>(n,m)P+", 0<n<m. 

Similarly, all bounded solutions of Un+i — ^nVn = rn,n < —1 are given by 

y- = <^>{n,0)r]_+ J2 G-in,m + l)rm, n < 0, G 7^(Po-"), (88) 

m<-l 

where 

, , / <l>(r^,m)P^^ m < n < 0, 
G-{n,m) = < J-/ \ n-,, (89) 
^ ' ^ [ -$(n, m)P^'', n < m < 0. ^ ^ 

By the exponential dichotomies the Green's functions satisfy 

\\G±{n,m)\\ < Ke-"l"-'"l, n,m G Z±. (90) 
With u^H'z '■= {u'z, Htl) we may write flHT]) in block operator form as 

(I T) tt) ^ (7) ■ <-) 
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By the bordering lemma [21 Appendix] the block operator is Fredholm of the same 
index as £ and, using (18^ . it is a linear homeomorphism in i°°{]R'') x R. Since 
£^ is a closed subspace of i°°(JR!') it suffices to prove that the unique solution 
(l/s, h) of (EI]) in ^""(R'^) x R satisfies the estimate (ED in case E C- 

Take the inner product of the first equation of (|ST]) with w^. Then and 
the normalization ||w;s||^2 = 1 show h = (if a, '"a)- Therefore, by (|T8|) . 

<Ce||rs|U$^e""l"la;„<Ce^^^||rsL. (92) 

With this we have Z'z := — hw'z G 7^(>C) by (18^ . Below we will construct a 
special solution G ^""(R'^) of £ys = -2z such that for some constant C > 

ll^zL < C'IkslU- (93) 

By ([HID and the solution of is given by 

{yw„h) = {y'z + cu'z,h), c = 7-(us,yz), h={wz,rz). (94) 

From the exponential decay ((Tj), (fT8|) and Lemma [161 applied to the kernel 

^-ai\n\ + \m\) < gag-a|„-m-l| ^ obtain with C = CeC^iTo, 

Then ([92])- ([94]) yield the assertion 

IbsL < (1 + c'Oll^aL + ItIII^zL < c'(i + c'O^lkalU + bllkzL- 

It remains to construct with = and ([93]) . We determine 77+ G 7^(Pq'"'*) 
and 77- G TZ{Pq^) such that y„ = y^, n > with y,;]^ from ( 1861) . and = ?/~, n < 
with ?/~ from ([88]) . and such that the definitions coincide at n = 0. The last 
condition holds if and only if 

m< — 1 m>0 

By ([87]) . ( [89]) . ([85]) the first sum on the right is in © Yi and the second sum 
is in Y_ © Yi while the left-hand side is in Fq © © Since 2;^ G 7^(>C) holds, 
equation ([95]) has a solution and thus Aq G ©y_. We conclude from ([85]) that 
= Pq^'^^o and r]_ = -Pq""Ao are the unique solutions of ([95]). With ([90]) and 
Lemma [TB] we estimate for n > 

WVnWoo;;' < KlI^sL 5^a;-ie-"l"-™^ila;,„ < KKo\\z^\\^. 

m>0 

An analogous estimate holds for n < and this completes the proof. ■ 
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Proof: (Lemma I14h We use Lemma [18] and construct an approximate right 

inverse 5+ of D(x,g)G^s = -D(a;,g)Gs(0, 0, 0, 0). For any £ G /(s) we define the 
interval J{€) = {£-, ...,£,...,£+} where right and left neighbors are given by 

^ ^ f cx), if e = maxl{s), 

+ \ i + N^, otherwise, ^ ' 

-oo, _ if £ = min/(s), 

max{£+ : i G I{s)J <£} + !, otherwise. ^ ' 

Note that the sets J{i),i E I{s) define a partitioning of li. In the following we 
consider N > N = 2N^: + 1 which implies £ — > ^+ — ^ > A^* • During the proof 
A^* will be taken sufficiently large. Given an element (^^,7) G £°°(JR!^) x i°°{s), 
we decompose 

r^= ^m'^TL, where (lj(,))„ = | J' ^^1^''' 

is the characteristic function of J{i)- Let Bq denote the solution operator of fl9T]) . 
then we set 

=5o(/3'l Jwrs, 7^) for ^ G J(s), (98) 
and define 5+ as a blockwise inverse via 

B+{r^, 7) = = ( 5^ (99) 

leiis) 

Using Lemma [T7] with the settings rii = i_ — i, Ur = i+ — i we obtain a bound 

Ib^L + \hi\ < C* (||/3'lj(f)rsL + < C* (llrslloo + • (100) 
Let us abbreviate the weights from ( 1791) . 

Wn,^ = tUn(a, ^+), n G £ G I{s). 
Then equation fl99|) and fllOOp lead to the estimate 

hnW < II 5Z ^n-^ll ^ '^^ (Ikzlloo + hiloo) C^n/, 

For 2e-''^* < 1 the last sum is bounded by 1 + 4e " and (11 00 1) yields 

||5+(rs,7)ll = ||ys||oo + ||/i||oo < (2C* + l + 4e-")(||rs||oo+ hiloo). (101) 
In the next step we show for A^* sufficiently large, 

||(rs,7) - A=^,,)G°i?+(rs,7)|| < ^(llr^loo + hlU), (102) 
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which by f lllSp . ( llOip gives the desired estimate 

||(Z}(.,,)G°)-^|| <2(2C* + l + 4e-»). 

We introduce (^^,7) = (^^,7) — D(^x,g)G%y'z,h). From (p9|) and the variational 
equation in fl98l) we have 

mGl{s) pGl{s) m£/(s) 

^ ^ ^ ^ 3;^— p7 0) fxi^n—rm^) yn—m- 

m£/(s) p£/(s) 

For n E there exists a unique £ G /(s) such that n G J(^). Define the neighbor- 
hood U{£) of £ by = {£, £, with left neighbor £ = sup{m G /(s) : m < £} 
and right neighbor £ = inf{m G /(s) : m > £} (as usual let £ = —00 and £ = 00 
if the sets are empty). Using the Lipschitz constant of fx and f|T3|) . fllOOp we 
obtain 

ll^nll ^ -^a; ^ ^ ] ll?/n-mll ^ ] ||3^n-p||+ ^ ] ll?/n-mll ^ ] H^^n-pl 



me/(s)\W(f) t^P&I{s) 



p&I{s)\{l/} p(.H^s)\{l/} 

We show that the terms in {. . .} are of order 0(e~°^*) so that the contraction 
estimate (11021) holds for the first component if A^* is sufficiently large. A critical 
term on the right-hand side is 

oj^^l\\xn-t\\ < c^e-'^\^-^\+<^-'^-) < Cee'^(^-^'-) < Cee"'^^*, £_<n< £+. (103) 

The term a; is handled analogously. Further, 



i^pel{s) l>p£l{s) i<p&I{s) 



' e>peiis) £<pei(s) 

The remaining terms allow similar estimates since n always lies in an exponential 
decaying tail of the weights and of the shifted homoclinic orbits. Finally, we use 
(HD and ([98D, (^M) for £ G /(s), 

he - = I XI Y yn-m)\ 

<C^eC^*(||rd|oo + ||7l|oo)Xe"°'"'" Y 
37 



We estimate the remaining sum by using fllOSp 

neZ ijtm&I{s) n=l- m&I{s)\U{e) 

+ c(5^e-"(""^)+ 5^ e^"(^-")). 



n=£4 



The last two terms are 0{e "^*) since i+—i> N* and > A^*. Furthermore, 
we have for i_ < n < i. 



_^a(m+— n) 

m(iI{s)\U{C) £<mel{s) £>me/(s) 



< a(£_-m_) , Mm+-e+) <- 

i<m&I{s) i>mel{s) 



Our final estimate is 



J2 e""'"~"c^„,,- < E e-"(^-")+'^("-^-) + ^ e-" 

n=e- n=e- n=l+l 



(n-£)+a(n-£_) 



< ^ + 



1 - e-" 1 - e-("-'^) ■ 
The term Ylf^=e^ e"°l"~^la;^ | is estimated in a similar way. This finishes the proof 

of (uni). 

We take the same B_ = as an approximate left inverse in Lemma [13 
It is convenient to require < 2a < Qq < a instead of (ISUjl . Given {y^^h) e 
£°°(R'=) X we set 

(r/s,t) = B+D(,,,)G'°(ys,/i), 
then it is sufficient to show 

\\m - vAoo + ||t - h\U < 0(e-("-")'^*)(||ys||oo + \\h\\oo). (104) 
For a fixed I G /(s) we consider first the case 

2/n = (n < and n > hm = {m ^ i), (105) 

and prove the stronger estimate (recall = 2iV* + 1) 

\\v^-yA^ + \ti-h\+ sup |t™|e"(l"^-^l-^)<(9(e-"^*)(||ysL + |/i£|), (106) 
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where Un = Un{a, £-,£+), see ( !79|l . Then we consider the case 

yn = (n^i^), h = (107) 

and prove the estimate 

\\V^ - y^L + sup |t„|e'^l'"-'-l < (9(e-'^^*)||ys|U, (108) 

mew, 

where a;„, = u{a, Let us first show that the general case fll04p follows from 

mm and (dug), with J^{£) = {£_ + !,.. .,£+}, J^{£) = {£_} we decompose 

(yz, ^) = X] /iml{m}) + (Iji(m)yz, 0)] 

mG/(s) 

and define 

B+D(^.^^g)G^^ [(lj0(m)yz,^ml{m}) + ( 1 Jl {m)l/Z , 0)] . 

»7ie/(s) 

Similar to the proof of (11021) we combine the local estimates (11061) and (llOSp . 

Ilz/n - ?7n|| = 0(6""^*) ^ [a;,,(a, m_,m+) (II ljo(„)ys L + |/im|) 

m£l(s) 

+ u;„,(a,m_,m_)||l{„_.}?/sL] 
= C(e"''^*)(||yz||oo + ||/i||oo) ^ [w„(a,m_,m+) + w„(a,m_,m_)] . 

= 0(e-'^^*)(|bdloo+||/^||oo). 

Using the exponential weights of in (I106P and (11081) leads to the same estimate 
for \\h — t\\oo and (I104p is proved. 

For the proof of ( 1106^ let ■jm = Vz) and note that (llOSp implies r„ = 

for n ^ J{£)- Hence r]^ = f3^yw, satisfies 

Vi+l - fxiXn, OWn + hiWn = T^+e + [fxiPn+lis], 0) - /x.(x„, 0)] ?7^. 

Using Lemma [T7| and a Lipschitz estimate for fx we can compare with the solution 
iy'^,t,) = 5o(/3^1j(,)rs,7^) of (EHD. This leads to 

Wr'yi - y^L + \te - hi\ = 0(e-"^*)(||2/^L + 

For the solutions {y^,tm) = 5o(/3™lj(m)'"z, 7m) = 5o(0,7m) with £ Lemma 
[T71 gives 

^+ 

||/3-"^2/^L + K\ < C*|7m| < CelbdU Yl = 

n=£_+l 
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Collecting the last two estimates we find that {rj'^^t) = (X]mG/(s) P~^y%^ {im)m£i{s)) 
(cf. (p9|) ) satisfies inequality f ll06p . 
In case condition f ll07p holds, let 

r„ = {D^F{p^{s),Q)y^)n {n E Z), 7™ = {m G /(s)), 

and note r„ = for n 7^ — 1 as well as 

re_ = -fxiP£_is),0)ye_, ri__i=yi_. (109) 

Moreover, 

= \uL-^yiA < Cee""l^-""^l||y,J|, m G /(s). (110) 

As in dM]) let (y^, t„) = 5o(/3'"lj(^)rs, 7m) for m G /(sj. Recall £_-l = 1+ from 
( 197|) where £ G /(s) is the left neighbor of £. For m ^ l. iwe have = /3~™7mMz, 
= by ([94D and therefore by (fTTOj) . 

Ilz/s L < C2|||/,_||e-"l^-'"lsupe-"l"-'-l+'^l""^-l = 0(e-("-'^)l^-'"l)||y,_||. (Ill) 

We introduce the weights = a;„(2a, £_ — — i). For m = i we use (11 101) 
and Lemma [T71 with w*, 

lb^L-<C*(ll/3'liWrdU* + l7.|||«dL0 .^^2) 
<C*(l + Ce-(°-2'^)l^-"^l) <C||y,„||. 

In a similar manner, 

\\yi\W < C\\yeJ (113) 
holds for the weights = Ci;„(2a, — £,£+ — £). The t- values satisfy 

It^l = \uj _J.(p^Js),0)y,J <Ce-"l^-^l||2/,_||, 

(114) 



In particular, this proves the t-estimates in f llOSp . Next we estimate the difference 
d'z = P~^y%+ P~^y%~y'E by using the exponential dichotomy on Z for the constant 
coefficient operator ^oZ/a = iyn+i — Ayn)n£'z, A = /x(0,0). From (I109p and the 
definition of y^,y^ we find 



dn+i - Adn = {fx{xn-e, 0) - A)y'^_^ + {fx{x^_i, 0) - A)y'^_ 

- teWn-e - tiw^_i +{A- f^ipe^ (s), 0))5„,^_y^_ = ^ T^. 



-i 

5 



j=l 



40 



For every term we show ||Ti|| = 0{e "'^*)cOn\\ye-\\ with weights Un = e ^^L In 
fact, our previous estimates (11121) - (11141) yield 

||Ti||u;-i<C||x„_,||e-'^l"-^-l||y^||.. 

< C||y,_||e-°l""^l"«l"-^-l < Ce-'^l^-^-l||y,J| = 0(e""^*)||y,_||, 

^sIK' < Clt.le'^l"-^-!-"!"-^! < Ce-("-'^)l^-^-l||y,J| = 0(e-'^^*)||y,J|, 
^4|K' < C|t,|e'^l"-^-l-"l"-^l < Ce-("-'^)l^"-^+l||y,J| = 0(e-'^^*)||y,J|, 
^5||a;„-'<C|| 5^ x,__„,||||y,_||=0(e-"^*)||y,_||. 

mG/(s) 

Since the operator Cq has a Green's function with an exponentially decaying 
kernel we infer from Lemma [16] 

Combining this with (II lip gives 



Jn—m I 



II 5^ C-n^-?/n|| < |Mn|| + 

< C(a;„e-"^* + e"'^!^-™!) |b,J| = 0(e-'^^*)a;„||y,J 

This finally proves the ^/-estimates in (11081) . ■ 

A Auxiliary results 

Lemma 18 (Banach Lemma) Let X,Y he Banach spaces and let A G F], 
B_,B^ e L[y, X] be bounded linear operators such that 

||/y - A5+II < 1, ||/y - fi-A|| < 1. 

Then A is a homeomorphism with 

Proof: Note that y = [ly — AB^)y + r has a unique solution y for every r G F. 
Then y satisfies 

- l-||Jy-A5+|| 
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and X = B^y solves Ax = r. To prove uniqueness, note that any solution x of 
Ax = r solves x = {Ix — B_A)x + B^r. Since Ix — B^A is also contractive the 
solution is unique and the estimates follow. ■ 
A key tool in the proofs of Lemma [I3] and Theorem is the following quan- 
titative version of the Lipschitz inverse mapping theorem, cf . [13] . 



Theorem 19 Assume Y and Z are Banach spaces, F G C^{Y, Z) and -F'(yo) 
for Hq E Y a homeomorphism. Let k, a, 6 > be three constants, such that the 
following estimates hold: 

\\F'{y)-F'{y^)\\ < < a < / Vy G ^^(yo), (116) 

\\F m) II 

||i^(yo)|| < (117) 
Then F has a unique zero y G Bs{yo) and the following inequalities are satisfied 

\\F\yr'\\ < yyeBsiyo), (118) 

a — K 

II2/1-I/2II < ^— ||F(yi)-F(y2)|| Vyi, 2/2 e 5^(1/0). (119) 
a — K 

We collect some well known results on exponential dichotomies from [22] . 
Denote by $ the solution operator of the linear difference equation 

yn+i = Anyn, « G (120) 

which is defined as 

{An-i . . . A^, for n> m, 
J, for n = m, (121) 

. . . forn<m. 

Definition 20 The linear difference equation (11201) with invertible matrices An G 
R*''''^ has an exponential dichotomy with data {K, a, P^, P^) on an interval 
J alt, if there exist two families of projectors P^ and P^ = I — P^ and constants 
K, a > 0, such that the following statements hold: 

Pn^{n, m) = ^{n, m)P^ Vn, m E J, 

II $(n,m)P4 II < Ke-"^"^ 

\/n > m, n,m E J. 

||$(m,n)P^|| < /s:e-"("-™) 

Theorem 21 (Roughness Theorem, cf. Proposition 2.10]) Assume that the 
difference equation 

yn+i = Anyn, An E R''''' invcrtibk, \\A~^\\ < M \fn E J 
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with an interval J C has an exponential dichotomy with data {K,a, P^, P^) . 
Suppose Bn G R^''^ satisfies ||-Bn|| < i> for all n ^ J with a sufficiently small b. 
Then An, + -B„ is invertible and the perturbed difference equation 

yn+l = {An + Bn)yn 

has an exponential dichotomy on J. 
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